Stability and quasinormal modes of the massive scalar field around Kerr black holes 
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We find quasinormal spectrum of the massive scalar field in the background of the Kerr black 
holes. We show that all found modes are damped under the quasinormal modes boundary conditions 
when /j,M is not large, thereby implying stability of the massive scalar field. This complements the 
region of stability determined by the Beyer inequality for large masses of the field. We show that, 
similar to the case of a non-rotating black holes, the massive term of the scalar field does not 
contribute in the regime of high damping. Thereby, the high damping asymptotic should be the 
same as for the massless scalar field. 
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I. INTRODUCTION 

Scalar field in a black hole backgound is a subject of 
active investigation QJ, yet the massive scalar field has 
been investigated onlyin very few studies as to its quasi- 
normal spectrum |2[ , Q , late-time behavior Q and scat- 
taring properties |fj| . Behavior of the massive scalar field 
in a black hole background is quite different from that 
of the massless one in many aspects: first, it demon- 
strates the so-called superradiant instability @, which 
is absent for massless fields. This super-radiant regime 
happens provided the real oscillation frequency satisfies 
the inequality |i?ecj| < ma/ (2Mr+), where m is the az- 
imutal number of the mode, and M is the mass of the 
black hole. Another research, indicating the instability 
of massive scalars was that by Detweiler Q where slowly 
growing instability was suggested for values [iM « 1, 
i.e. in the regime, when particle Compton wavelength is 
much larger then the size of the black hole (/i is the mass 
of the scalar field). Second, it may cause existence of 
infinitely long- living modes called quasi-resonances @. 
Finally, at asymptotically late times the massive fields 
show universal behavior independent on spin of the field 
0. Note, that scalar field with the mass term can be 
also interpreted as a self-interacting scalar field within 
regime of small perturbations fTfH . One more recent ap- 
plication of dynamic of the massive scalar field in a black 
hole backgound comes from different brane-world theo- 
ries where massless fields can be considered as effectively 
massive fields on the brane. 

Yet, the strongest motivation to make an extensive 
search for quasinormal modes of massive scalar field is 
that until now the stability for massive scalar field pertur- 
bations in the Kerr background is not studied under non- 
super-radiant (quasinormal) boundary conditions. Usu- 
ally, stability of spherically symmetric black holes, like 
Schwarzschild or Reissner-Nordrstrom ones, can be rela- 
tively easily proved by using positivity of the conserved 
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energy. On the contrary to the Schwarzschild case, Kerr 
black holes have negative energy density inside the ergo- 
sphere. Therefore, fields can grow in parts of the space- 
time, leading to instability. The only result on stability 
of massive scalar field in Kerr space-time, we are aware 
of, is that of Beyer [llj], where it was proved that the 
modes of field with the mass obeying the inequality 

\m\a I 2M 1$ 
2Mr + y r + r+ 

are stable. Note, that the inequality found by Beyer is 
not compatible with conditions for super-radiance, so it 
does not contradict super-radiant instability of massive 
fields. In addition, the above inequality does not mean 
that fields of lower mass are unstable, and Beyer sug- 
gested the ways to lower the inequality 11]. If such 
an instability exists one could find an accurate thresh- 
old value of field mass y, when the instability begins. 
This could be done by an extensive computing of the 
quasinormal modes and by finding of the growing ones. 
Some low laying WKB quasinormal modes were found 
by Simone and Will j3j, yet neither fundamental modes 
(£ = n = 0, I is multi-pole number), nor modes with 
n > i can be found with trustable accuracy by WKB 
approach O O E3, 

even despite very good accuracy 
of WKB method for low overtones [l^. Therefore, being 
limited by WKB accuracy, one cannot judge about stabil- 
ity and high damping asymptotic in this case. Recently 
the quasinormal damping and late time regimes in time 
domain were investigated by L. Burko and G. Khanna 
111 and no instability was observed in the considered 
range of values of field mass fi ||. 

Because of the above reasons, in this paper we used the 
accurate Frobenius method to make an extensive search 
for quasinormal modes of the massive scalar field around 
Kerr black holes. The paper organized as follows: Sec. II 
introduces the main formulas for the wave equations and 
the Frobenius procedure. Sec. Ill gives numerical data 
for found quasinormal modes, and Sec. IV discusses the 
obtained results and the possibility for (un)stability. 
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II. BASIC EQUATIONS 

The background metric in the Boyer-Lindquist coordi- 
nates t, r, 0, (f> has the form: 



AaMr sin 2 6 



1 dt 2 H dtdch - — J ~ 2 

S / S A 



,io 9 2Ma 2 r sin 2 i 
2 ' r + a H 1 ™ a ^- k 



dtd(f> — —dr 

sin 2 6d<j) 2 , (2) 



where, M is the mass of the black hole, 
A = r 2 +a 2 - 2Mr, 

and 

£ = r 2 + a 2 cos 2 0. 

The massive scalar field obeys the equation: 

□$ + /i 2 $ = 0. (3) 

The radial wave equation ($ = e~ luJt e lm4, R(r)S((p, 9)) is 
well-known: 



d_ ( A dR(r) 



^--A-|iV )/?(./ )= 0. (4) 



where 

K = uj(t 2 + a 2 ) - am. (5) 
The separation constant A satisfies the equation 



1 d 



o dS(0) \ 
smOde V d9 ) 



sin 



(6) 



-a z u z sin 2 6» - a> 2 cos 2 (9 + 2amu + Xj S(9) = 0, 

and can be found numerically for any given u> |16| . 

Since A(w) can be obtained numerically, we are able to 
solve the equation for the radial part @, using the con- 
tinued fraction method |17| . In order to satisfy the quasi- 
normal mode boundary conditions, which corresponds to 
purely in-going waves at the event horizon and purely 
outgoing waves at spatial infinity, we require that 



R(r — > oo ) oc exp(zxr)r 1. 



R(r — » r+) oc (r - r + ) a = 



= (a; 2 + x ) 
a->r + (r^_ + a 2 ) — mar + 



r\ — a 2 



where r + marks the event horizon: 

r+ = M + VM 2 - a 2 , 

and 



X = ±\/uj 2 - /i 2 



The sign for x should be chosen in order to remain in the 
same complex quadrant as to. 

The equation has three singularity points: the spa- 
tial infinity (r — oo), the event horizon (r = Tj r ) and the 
internal horizon (r = r_ = a 2 /r + < r + ). Thefore the 
appropriate Frobenius series has the form: 

. / r a 2 \ l<r_ / rr + - r 2 + x 
E(r) = exp(jxr) r ^ I 



r7' + — a 2 



•(7) 



And the coefficients satisfy the three-term recurrence 
relation: 

a n a n+ i + f3 n a n + ~f n a„~i = 0, n > 0, 70 = 0, (8) 

where a n , f3 n , j n can be found in analytical form. We do 
not present them here because they have rather comber- 
some form. 

By comparing the ratio of the series coefficients 



In 



Otn-l 



Ot n (3n- 



a 



n 



Pn 



an-2ln-l 

1 /3„_ 2 — a„_37»i-2/. 
7n+l 

■ 1 +17ti + 2 



Pn 



(9) 



+ 1 P,i + 2-a n + 27n + 3/- 



we obtain an equation with a convergent infinite contin- 
ued fraction on its right side: 



Pn 



Pn 



a„-27n-l 
/3n-2 — ffln-STn-a/— 

a n 7n+i 



(10) 



/3„ + 2 -Q„ + 27n + 3/- 



that can be solved numerically by minimizing the ab- 
solute value of the difference between its left and right 
sides. The equation (|10|l has infinite number of roots, but 
the most stable root depends on n. Generally the larger 
number n corresponds to the larger imaginary part of the 
root w. 

Note, that the case under consideration allows to use 
the Nollert procedure |18| . in order to improve conver- 
gence of the infinite continued fraction, that is useful for 
searching roots with a very large imaginary part. 



III. NUMERICAL RESULTS 

Since we are limited by Kaluza-Klein equation, we do 
not take into consideration the effect of the mass of the 
scalar field onto the black hole. Thereby, we can consider 
only small values of mass of the field /1, in comparison 
with the mass of the black hole M, i.e. /jM << 1. Yet, 
as large values of /i are possible in brane- world theories, 
we shall touch this case as well. 

The quasinormal modes of the Kerr black hole depend 
on several parameters: the multi-pole number £, the over- 
tone number n, the azumutal number m, and on the black 
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Table I: Values of the quasinormal frequencies for the fundamental mode n = 0, £ — m — for different values of mass fx, and 
rotation a. 



l = m = Q n = 0.1 ll = 0.2 /i = 0.3 



a 


Re(coo) 


-Im(tJo) 


Re(aj ) 


-Im(wo) 


Re(cj ) 


-Im(tJo) 





0.112361 


0.096823 


0.116207 


0.075358 


0.122060 


0.043605 


0.1 


0.112447 


0.096741 


0.116311 


0.075312 


0.122442 


0.043495 


0.2 


0.112707 


0.096490 


0.116625 


0.075167 


0.123554 


0.043800 


0.3 


0.113140 


0.096050 


0.117153 


0.074910 


0.123993 


0.044905 


0.4 


0.113747 


0.095387 


0.117903 


0.074511 


0.123805 


0.044021 


0.5 


0.114521 


0.094442 


0.118883 


0.073988 


0.126101 


0.045025 


0.6 


0.115442 


0.093120 


0.120093 


0.073056 


0.128069 


0.044110 


0.7 


0.116442 


0.091261 


0.121514 


0.071782 


0.130044 


0.044193 


0.8 


0.117301 


0.088592 


0.123033 


0.069849 


0.129843 


0.043152 


0.9 


0.117181 


0.084742 


0.124136 


0.066839 


0.135784 


0.041884 


0.95 


0.115848 


0.082697 


0.123910 


0.064981 


0.136616 


0.042938 


0.99 


0.114478 


0.082285 


0.123197 


0.064010 


0.134994 


0.040595 


0.995 


0.114388 


0.082240 


0.123144 


0.063928 


0.136560 


0.040180 



hole and scalar field parameters: M, a, /i. Note that m 
can take values —£, — £ + 1, ...,£ — 1, £. From here and on 
we measured /u, and a in units of black hole mass M. 

Let us look at tables I-III: the damping rate determined 
by the imaginary part of the uj is decreasing considerably 
with increasing of mass /i. While the real oscillation fre- 
quency is slightly increased with grwoingof \x. This hap- 
pens also for Schwarzschild black hole 2]. Note, that it 
is the fundamental mode n = that contributes mainly 
into the ringing signal. 



To check the correctness of the applied procedure we 
considered particular cases of massless scalar field in 
the Kerr background and of the massive field in the 
Schwarzschild background. For these two cases the quasi- 
normal modes are well-known for example from the pa- 
pers 0], 0, , One can see, for instance in the 
table III that QNMs in case of /u, = 0, coincide with cor- 
responding data in of 0, 0, 0- Also data in tables 
I- IV for a = coincide with modes found in 



As the number of the overtone is increased, the QN 
frequency of the massive field becomes closer to that of 
the massless field. For the case of Schwarzschild black 
hole, we proved in |8( , both analytically and numerically, 
that this tendency of QN modes to diminsh difference 
stipulated by mass fj, 7 ends with a universal high damp- 
ing asymptotic, which is independent on the mass of the 
field. Thus, the massive term /j, does not contribute to 
high overtone behavior of Schwarzschild black holes. We 
shall show numerically that the same is true for Kerr 
black holes. Indeed, from Table II, one can see that al- 
ready for relatively small overtone number, the difference 
between quasinormal modes with different values of mass 
/i becomes negligible. This is important for us, because 
now, we are in position to say that it is enough to study 
low overtone behavior to judge about stability of Kerr 
case within our numerical approach. 

From Tables I-V and figures (1-2), we see that all found 
modes for different values of parameters I and m are 
damped in the regime of small /i. Larger fj, (see Fig. 
3) leads to considerable decreasing of the damping rate 
and finally to the same infinitely long living modes (i.e. 
modes with infinitesimaly small imaginary part) found 
in Q for Schwarzschild black holes and called quasi- 



Table II: 


QNMs for 


fi — 0.3 and ii 


= 0.1, a = 0.5, 


e = m = i. 


n 


M 


= 0.3 


M = 


0.1 


1 


0.330086 


- 0.267718 i 


0.325510 - 


0.288297 i 


2 


0.292386 


- 0.489967 i 


0.296315 - 


0.501634 i 


3 


0.266009 


- 0.720273 i 


0.270472 - 


0.726246 i 


4 


0.245565 


- 0.952427 i 


0.249182 - 


0.955699 i 


5 


0.228436 


- 1.184966 i 


0.231248 - 


1.186916 i 


10 


0.164599 


- 2.343964 i 


0.165873 - 


2.344194 i 


15 


0.117999 


- 3.499102 i 


0.118599 - 


3.499165 i 



resonances. Therefore, when increasing /i, the strict line 
in Fig. 3 does not have continuation into negative half- 
plane. We see here that, the larger mass of the field 
fi is, the more low laying modes will change into quasi- 
resonances. Meanwhile, the rest of the spectra consists 
of ordinary modes (i.e. modes with finite time of damp- 
ing). Moreover, at sufficiently high overtone number, 
the difference between massive and massless case will be 
neglected (see Table II). This behavior mimics the QN 
behavior of non-rotating black holes. Yet, potential for 
massive scalar field around non-rotating black holes is 
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-0.0953 
-0.09535 

-0.0954 
-0.09545 

-i.0955 



Figure 1: Imui as a function of azimutal number m for 
a = 0.3, n = 0, fj, = 0.1. 



Figure 2: Jmoj as a function of azimutal number m for 
a = 0.3, n = 0, fj, = 0.1. 



definitely positive everywhere outside black hole, while 
for Kerr solution, the potential depends upon complex 



frequency and such a stightforward analysis is impossi- 
ble. 

Now we are in position to state that the high damping 
asymptotic for quasinormal modes of massive scalar field 
is the same as for massless one, i.e. the real part increases 
with m according to the numerical law described in |2Cj . 

To exclude possible instability we should search for 
QNMs for different values of £ and to. As one can see 
from representative plots Fig. 1, at large £ the imagi- 
nary part of the QN frequencies remains bounded and, 
thereby, does not show any tendency to instability. The 
dependence of the real osicillation frequency on the az- 
imutal number to can be found on Fig. 3. 

After careful investigating of region of small values 
/xM, we conlcude that there are no unstable modes. As 
to the case of intermediate /aM, our extensive search of 
quasinormal modes implies that there are no unstable 
modes at least as far as fiM is not much larger then 
1. We certainly cannot find quasinormal modes for any 
large value of fi. This happens because the numerical 
procedure used here requires slow changing of [a, when 
searching for new modes, in order not to jump occassion- 
al into another overtone. Thereby, considerably large 
values of \i can be reached only for very long time of 
computation. 
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Table III: Values of the quasinormal frequencies for the fundamental mode n = 0, I = 1, m 
and rotation a. 



for different values of mass fj,, 



l:m = 



/Lt = 0.1 



H = 0.2 



0.3 



a 


Re(tJo) 


-Im(wo) 


Re(uj ) 


-Im(wo) 


Re(ioo) 


-Im(wo) 





0.297416 


0.094957 


0.310957 


0.086593 


0.333777 


0.071658 


0.1 


0.297602 


0.094884 


0.311128 


0.086544 


0.333924 


0.071646 


0.2 


0.298164 


0.094661 


0.311646 


0.086392 


0.334366 


0.071608 


0.3 


0.299116 


0.094273 


0.312523 


0.086124 


0.335116 


0.044905 


0.4 


0.300478 


0.093692 


0.313780 


0.085717 


0.336190 


0.071408 


0.5 


0.302285 


0.092873 


0.315446 


0.085132 


0.337622 


0.711202 


0.6 


0.304579 


0.091742 


0.317564 


0.084308 


0.339444 


0.070869 


0.7 


0.307414 


0.090182 


0.320185 


0.083141 


0.341707 


0.070335 


0.8 


0.310836 


0.087989 


0.323357 


0.081461 


0.344465 


0.069472 


0.9 


0.314815 


0.084810 


0.327070 


0.078967 


0.347735 


0.068059 


0.95 


0.316919 


0.082696 


0.329061 


0.077284 


0.347735 


0.068058 


0.99 


0.318576 


0.080718 


0.351004 


0.066054 


0.351004 


0.066054 


0.995 


0.318779 


0.080453 


0.351189 


0.065918 


0.351189 


0.065918 



IV. CONCLUSION 



Within numerical approach we managed to compute 



In the present paper we have considered the mas- 
sive scalar field perturbations for the Kerr black holes. 
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Table IV: Values of the quasinormal frequencies for the fundamental mode n = 0, £ = 1, m = 1 for different values of mass /x, 
and rotation a. 

£=l;m = l M = /x = 0.1 p = 0.2 



a 


Re(tJo) 


-Im(wo) 


Re(a>o) 


-Im(o;o) 


Re(coo) 


-Im(wo) 





0.292936 


0.097660 


0.297416 


0.094957 


0.310957 


0.086593 


0.1 


0.301045 


0.097547 


0.305329 


0.095029 


0.318274 


0.087228 


0.2 


0.310043 


0.097245 


0.314119 


0.094920 


0.326433 


0.087709 


0.3 


0.320126 


0.096691 


0.323981 


0.094569 


0.335621 


0.087979 


0.4 


0.331567 


0.095792 


0.335181 


0.093883 


0.346095 


0.087950 


0.5 


0.344753 


0.094395 


0.348105 


0.092714 


0.358230 


0.087478 


0.6 


0.360285 


0.092243 


0.363345 


0.090805 


0.372594 


0.086320 


0.7 


0.379159 


0.088848 


0.381888 


0.087678 


0.390141 


0.084014 


0.8 


0.403273 


0.083132 


0.405606 


0.082262 


0.412675 


0.079526 


0.9 


0.437234 


0.071848 


0.439045 


0.071342 


0.444549 


0.069737 


0.95 


0.462261 


0.060091 


0.463691 


0.059825 


0.468050 


0.058968 


0.99 


0.493423 


0.036712 


0.494284 


0.036756 


0.496939 


0.036879 


0.995 


0.490424 


0.073874 


0.490872 


0.073846 


0.492223 


0.073740 



Table V: Values of the quasinormal frequencies for higher overtones for fiM = 0.1, I = 1, m = 1. 



i=l;m = l 




n= 1 




n = 2 




n = 3 


a 


Re(wo) 


-Im(wo) 


Re(cj ) 


-Im(wo) 


Re(wo) 


-Im(wo) 





0.264689 


0.302851 


0.228650 


0.538599 


0.202558 


0.787632 


0.1 


0.274260 


0.301711 


0.239156 


0.534670 


0.212898 


0.780859 


0.2 


0.284872 


0.299936 


0.250883 


0.529466 


0.224537 


0.772046 


0.3 


0.296738 


0.297332 


0.264079 


0.522633 


0.237749 


0.760644 


0.4 


0.310140 


0.293608 


0.279072 


0.513640 


0.252894 


0.745820 


0.5 


0.325510 


0.288297 


0.296315 


0.501635 


0.270471 


0.726246 


0.6 


0.343441 


0.280603 


0.316481 


0.485169 


0.291244 


0.699664 


0.7 


0.364932 


0.269019 


0.340644 


0.461505 


0.316511 


0.661835 


0.8 


0.391790 


0.250187 


0.370780 


0.424542 


0.348980 


0.603557 


0.9 


0.428112 


0.213833 


0.411782 


0.355839 


0.382103 


0.644425 


0.95 


0.441634 


0.287111 


0.431664 


0.399149 


0.382103 


0.644425 


0.99 


0.473604 


0.275655 


0.468124 


0.395139 


0.460787 


0.576199 


0.995 


0.480613 


0.293136 


0.477664 


0.381687 


0.471122 


0.604783 



-lm(«) 




»=!!«(") 



0.08 0.085 0.09 0.095 0.1 0.105 0.11 



Figure 3: Reui vs Imui for n = 1, a = 0.5, I = m = for 
different values of n 



quasinormal modes for values \xM not much larger then 
1 . In this regime all found low laying modes are damped. 



We have shown that similar to the Schwrazschild case, 
the massive term does not contribute into the high damp- 
ing part of the spectrum. Therefore we can conclude that 
for not large values of /iM there is no instability for mas- 
sive scalar field under the quasinormal mode boundary 
conditions. From the inequality (1) it follows that for a 
given black hole mass M and rotation a, for any large 
there is sufficiently large m, so that inequality is not 
performed. Therefore we sould have instability for suffi- 
ciently large values of m, even for large /x. Yet, we have 
whon in Fig. 1, that QN spectrum for larger values of m 
does not show any tendency to instability, because the 
Irnto for increasing values of m remains bounded within 
some region of negative values corresponding to damp- 
ing. Thefore inequlity (1) found in |ll| means stability 
for fields with sufficiently large values of fi. 

In addition, we have shown that the qusi-resonances, 
which exist for the Schwrazschild metric, exist also for 
Kerr black holes. 
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